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Abstract 

A newly developed global dynamic model was 
used to simulate the dynamics of a gear noise rig 
at the NASA Lewis Research Center. Experi- 
mental results from the test rig were used to 
verify the analytical model. In this model, the 
number of degrees of freedom of the system are 
reduced by transforming the system equations of 
motion into modal coordinates. The vibration of 
the individual gear-shaft systems is coupled 
through the gear-mesh forces. A three- 
dimensional bearing model was used to couple the 
casing structural vibration to the gear rotor 
dynamics. The system of modal equations is 
solved to predict the resulting vibration at sever- 
al locations on the test rig. Experimental vibra- 
tion data were measured at several running 
speeds and were compared to the predictions of 
the global dynamic model. There was excellent 
agreement between the analytical and experimen- 
tal vibration results. 

Nomenclature 

{A} rotor modal displacement of (X, 0 X ) 
{A c } casing modal displacement of (X c , X c ^) 
{B} rotor modal displacement of (Y, 0 y ) 


W 

rotor modal displacement of (Y c , y c ^) 

[C b ] 

bearing damping matrix 

py 

[*I T icy [*] 


W T [cj s»d 

pj 

casing structure damping matrix 

Pci 

l*J T [CJ f»J 

{D} 

rotor modal displacement of Z 

{D c } 

casing modal displacement of (Z c , Z Q g) 

{D t } 

rotor modal displacement of 

F 

gear force 


external and mass imbalance 
excitations 

(F c (t)} 

force acting on casing structure 

(F G (t)} 

nonlinear gear forces, through gear 
mesh coupling 

{F g (t)} 

shaft bow force = [KJ {W r } 

Pa) 

rotor angular acceleration 


1 


Pa) 

[o T | 

[G t I 

P*1 

w 

[KJ 

[KJ 

[Sj 

[KJ 

[KJ 

(Ml 

[MJ 

R 


l»| T [G a ] [*] 

gyroscopic forces 

l*l T |Gvl (*] 

1/2 {[KJ + (K by ]> 

l*] T [K a 1 1*1 

bearing axial and lateral cross-coupling 
stiffness 

l*] T [KJ [*1 
l*J T (KJ (*«] 

casing structure stiffness 
shaft bow stiffness 
mass matrix of rotor 
mass matrix of casing structure 
pitch radius of gear 


{[W]> 


generalized displacement vector of 
rotor 


{[W c ]} generalized displacement vector of 
casing 

{W r } genearlized displacement vector of 
rotor 


X,Y,Z displacement vectors 


[*l 

W 

M 


friction coefficient between the gear 
tooth surface 

orthonormal mode 

orthonormal mode of casing 

critical speeds of each rotor 

critical speeds of casing 


Introduction 


Large vibrations in gear transmission systems 
causes excessive wear and crack formation in gear 
teeth, which results in premature gear failure. 
With the need for higher operating speeds and 
power from transmission systems, the problem of 
excessive vibration becomes even more critical. 

In order to insure smooth and safe operation, it is 
necessary to understand the dynamics of the gear 
transmission system. 

Two areas of research in the dynamics of gear 
transmission systems are (1) analytical simula- 
tion and (2) experimental testing. There is a 
great deal of literature on the vibration analysis 
of a single gear stage. 1 ’ 6 Some work has been 
done on multistage gear vibration, 7 ’ 9 but very 
limited work 10 " 11 has been done on the dynamic 
analysis of gearbox vibration. Considerable effort 
has been devoted to experimentally studying gear 
dyanmics 12 ’ 14 and localized vibration effects on 
gear teeth. 15 A few studies have been conducted 
to correlate analytically predicted and experi- 
mentally measured gearbox vibrations. 


a angle of orientation 

0 Eq. (5) 

$ generalized displacement 

torsional rotational vector 

9 , 9 lateral rotational vectors 

x y 


This paper correlates the experimental results 
obtained from the test rig at the NASA Lewis 
Research Center with predictions from an analyt- 
ical model developed by using the modal synthe- 
sis method. 7 The major excitations of the rotor 
system include mass imbalance, shaft residual 
bow, nonlinear gear-mesh forces 16 and gear-tooth 
frictional effects. 14 The vibratory motion 
between the rotor and the casing is coupled 


2 


through the support bearing in the lateral and 
axial directions. 1 Gearbox mode shapes and 
vibration predictions from the analytical model 
are compared to those obtained from experimen- 
tal testing. 


Analytical Procedure 

Development of Equations of Motion 

The equations of motion for each ^ear-shaft 
system can be written in matrix form 7 as 

[M] (w} + (G v l {w} + [G a ] {W} + [C b ] (W - W c } 
+ [K b ]{w-w c } + [iy(w- w r } 
= {F(t)} + {F G (t)J 

( 1 ) 


where [M] is the mass matrix of the rotor (iner- 
tia), {w} is the generalized displacement vector 
consisting of the three displacement vectors 
X,Y,Z with the corresponding lateral 0 X , d y , and 
torsional 0 t rotational vectors as 


{w} 


(X) 

(*z) 

(Y) . 

(Z) 

(•*) 


( 2 ) 


and where [G ] is the gyroscopic force; [G A ] is the 
rotor angular acceleration; [C b ] is the bearing 
direct and cross-coupling damping, [K b ] is the 
bearing axial and lateral cross-coupling stiff- 
ness; 11 {W — W c } is the casing vibration; [Kj is 
the shaft bow stiffness; {W — W r } is the shaft 
residual bow; (F(t)} is the external and mass 
imbalance excitation; and {F G (t)} is the nonlin- 
ear gear force through gear-mesh coupling. 


For a multiple gear-shaft system, the equa- 
tions of motion (1) are repeated for each shaft. 


The motions of the individual shafts are coupled 
to each other through the gear-mesh forces 
(Fig. 1) and to the casing through bearing stiff- 
ness [K b ] and damping [C b ]. 

The equations of motion for the casing can be 
written as 

[Mj{W c } + [C b ]{w c - W} + [K b ||w c - W} 

+ [Cj{wJ + |KJ{W.} = (F c (t)} 

(*) 

where [M ] is the mass matrix of the casing struc- 
ture; {W c } represents the generalized displace- 
ment vector of the casing structure: 


w = ■ 


[Cj is the casing structure damping matrix; [KJ 
is the casing structure stiffness; and (F c (t)} is 
the force acting on the casing structure. The 
nonlinear gear forces for the k th individual ^ear- 
shaft system, using nonlinear gear stiffness 10 and 
gear-tooth friction, 14 are given as 
x-force 

Fcxk = 52 X tki [“^ci^ci ~ ^ck^ck 

i=l, i *k 

+ ( X ci ” X ck) cos a ki (5) 

+ (Y ci - Y ck )sin a ki ] 

[cos ay + (sign) (p) (sin a ki )] 

y-force 


(X c ) 

(Xci) 

(Y c ) 

(Yci) 

CZe) 

(Zci) 
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^Gyk — E X tki [ ^ci^ci ^ck^ck 
i=l,ii*k 

+ ( X ci - X ck) cos a ki (6) 

+ (Y ci - Y ck )sin a u ] 

[sin a ki + (sign) (ft) (cos a ki )] 
and torsional 

^Gtk == 5^ ^ck ( X tki [( — ^ci^ci — ^ck^ck) 

i=l,i #k 

+ ( X ci ■- X ck) cos a ki 

+ (Y c i r Y ck )sin Qjy]} 

(7) 


where R is the pitch radius of the gear and fx 
is the coefficient of friction between the gear 
tooth surface and “SIGN” is the unity sign func- 
tion to provide the sign change when the mating 
teeth pass the pitch point. 14 


Modal Transformation 

The equation of motion for the undamped 
rotor system is 

[M|{WH-[[K,1 + [K a |]{W} = 0 (8) 

with the average bearing support stiffness from 
the x-y direction 7 given as 

[KaI = -‘{[K b ,l + [K by ]} (9) 

The orthogonality condition for the orthonormal 
modes [$] are 

[*] t [M][$] = [I] 

and 

[*] T [K 8 + K a ][$] = [c 2 ] 


( 10 ) 

( 11 ) 


where a; is the rotor critical speed. Similarly, a 
set of orthogonality conditions can be derived for 
the casing equation of 

[MJ {W c } + [CJ {W c > + [KJ {W c > = 0 (12) 

with the orthonormal mode [$ c ] such that 

W T KlW = (i] < 13 > 


W T [c c ]W = PJ 


M T ik c ] f»j = < i5 > 

where u) z is the casing critical speed. Using 
modal transformation 17 by letting 


{W} 


W (A) 
[**»}{ A) 

. [^y] (»} . 
l*j 9 } {B} 

M {D} 

w m 


(16) 


and 


[*«] (A c } 


PaJ{Ae} 


{W c } 


J[^cy] <B C } I 


[^cytfl (B c } 


(17) 


f*cJ (D c > 


l*cJ {D C }J 


the equations of motion for the rotor (Eq. (1)) 
can be transformed as 
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[I] {Z} + P,HZ> + 15*1 {z> + p b HZ} 

-(K b -K A ]{Z> + |*l T [C b lltJ{Z c } (18) 
+ |«*]{Z> - [*l T [K b ][*J{Z c > 

= 1*1 T {F(t) + F G (t) +F,(0) 


and for the casing (Eq. (3)) as 

[i|{zj + icj{zj + l^| {zj 

+ lK b ]{Z c > + p b l {ZJ - [#J T 
[K b l[*l<Z> - [*J T [C b ][*]{Z> 

= [*cl T Wt» 


where 


{Z} 


' ' 

{A} 

{B} 

A > 

{D} 

{DJ 


{Z c > 


{A c } 

‘ {B c } • 

< D c>. 


(19) 


( 20 ) 


Experimental Study 

The gear noise rig (Fig. 2) was used to mea- 
sure the vibration, dynamic load, and noise of a 
geared transmission. The rig features a simple 
gearbox (Fig. 3) containing a pair of parallel axis 
gears supported by rolling element bearings. A 
150-kW (200-hp), variable-speed electric motor 
powers the rig at one end, and an eddy-current 
dynamometer applies power-absorbing torque at 
the other end. The test gear parameters are 
given in Table 1. 

Two sets of experiments were performed on 
the gearbox; (1) experimental modal analysis and 
(2) dynamic vibration measurements during oper- 
ation. In the experimental modal analysis, modal 
parameters, such as system natural frequencies 
and their corresponding mode shapes, were ob- 
tained through transfer function measurements 
by using a two- channel, dynamic signal analyzer 
and modal analysis software. For this experi- 
ment, 116 nodes were selected on the gearbox 


housing. The measured three-dimensional mode 
shapes are presented in Fig. 4. The dynamic 
vibration measurements consist of data collected 
from accelerometers placed at three of the nodes 
on the surface of the casing. They were chosen 
such that vibration was measured in all three 
directions, X,Y, and Z. A dynamic signal analyz- 
er was used to compute the frequency spectra of 
the vibration. The experimental frequency spec- 
tra are shown in Figs. 5(a), 6(a), and 7(a) for the 
X-, Y-, and Z-directions, respectively, at a num- 
ber of different operating speeds. 


Discussion and Correlation of Results 

The measured mode shapes, shown in Fig. 4 
represent the major vibration modes of the gear 
noise rig in the 0- to 3-kHz region. Although 
these modes are only a small portion of the total 
modes of the system, they represent the major 
portion of the total global vibration of the sys- 
tem. In order to produce a compatible analytical 
simulation of the test apparatus, a similar set of 
modes was predicted by a finite-element model of 
the gearbox structure. This model serves as the 
basis for predicting casing vibrations in the over- 
all global dynamic model. Of the 25 modes found 
by the analytical model in the 0- to 3-kHz- 
frequency region, the 8 dominant modes were 
used to represent the gearbox dynamic character- 
istics. These simulated modes are shown in 
Fig. 8. The natural frequencies of the predicted 
modes are within 5 percent of the measured 
modes (Table 2). Also, the predicted mode 
shapes are very similar to the experimental 
modes shapes (Fig. 4). The good correlation of 
between the analytical model and the measure- 
ments confirms the accuracy of the dynamic 
representation of the test gearbox using only a 
limited number of modes. 

For the dynamic study of the gearbox vibra- 
tion, it was found that during a slow roll (low- 
speed run) of the gear-rotor assembly, a substan- 
tial residual bow (or eccentricity in the sleeve 
assembly) exists in the rotor system as shown by 
its large orbital motion in Fig. 9. Figure 9(a) 
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represents the orbit of the driver rotor at the 
gear location, and Fig. 9(b) represents that of the 
driven rotor. Note that the circular orbit in the 
driver rotor at low speed represents the residual 
bow deformation of the rotor. The elliptical orbit 
in the driven rotor is attributed to a combination 
of the residual bow effects and the vertical gear 
force from the torque of the driving rotor. In 
order to analytically simulate the influence of this 
effect, a residual bow of 2 mils (0.05 mm) is 
incorporated into the numerical model (Eq. (l)). 

The frequency spectra of the analytically pre- 
dicted casing vibration in the X-, Y-, and Z- 
directions are presented in Figs. 5(b), 6(b), and 
7(b), respectively. As seen in Fig. 5(a), the ex- 
perimental casing vibration in the X-direction 
shows a major vibration component at the gear- 
mesh frequency (28 times shaft speed) at each 
rotational speed. A closer examination of this 
excited frequency component shows that two 
major vibration peaks occur at the running 
speeds of 1500 rpm (at a tooth pass frequency of 
700 Hz) and 5500 rpm (at a tooth pass frequency 
of 2560 Hz). These peaks are a result of the 
tooth pass frequency exciting two of the major 
natural frequencies of the housing, namely the 
658- and 2536-Hz modes. The presence of other 
modes can be seen; however, the 658- and 
2536-Hz modes, when excited by the gear-mesh 
frequency, dominate the spectra. 

Comparing the predicted vibration spectra 
with the measured spectra reveals that, although 
the actual amplitude values did not always agree, 
the general trends of the spectra were very simi- 
lar. The predicted vibration spectra of the hous- 
ing in the X-direction is shown in Fig. 5(b). A 
comparison of Figs. 5(a) and (b) shows that the 
predicted amplitude at the gear-mesh frequency 
at 1500 rpm is only 3 percent above the measured 
value. The comparison at 5500 rpm is not that 
close; the predicted amplitude is 38 percent below 
the measured value. If trends are compared, the 
predicted spectra show the same gear-mesh, 
frequency-induced excitation of the 658- and 
2536-Hz modes as that found in the measured 


spectra at the running speeds of 1500 and 
5500 rpm, respectively. 

Figures 6 and 7 compare the predicted mea- 
sured housing vibration spectra in the Y- and Z- 
directions, respectively. The results of the com- 
parison are the same as those presented for the 
housing vibration in the X-direction (Figs. 5(a) 
and (b)). Acutal values of the components in the 
spectra were not always in good agreement; how- 
ever, the general trends for the predicted and 
measured housing vibration spectra were very 
similar. Also, as seen in Fig. 7(b) at the 
1500- rpm speed, the model predicts the second 
and third harmonics of the gear- mesh frequency. 
As shown in Fig. 7(a), the measured vibration 
confirms the presence of these two harmonics at 
the 1500-rpm running speed. 

Conclusions 

A newly developed global dynamic model was 
used to simulate the dynamics of a simple 
transmission system. Predicted casing vibrations 
were compared to measured results from the gear 
noise test rig at the NASA Lewis Research Cen- 
ter. The conclusions of this study are summa- 
rized as follows: 

L The dynamics of the housing can be accu- 
rately modeled with a limited amount of analyti- 
cally predicted, experimentally verified vibration 
modes of the structure. 

2. The global dynamic model is capable of 
including in the analysis the effects of shaft resid- 
ual bow or eccentricity. 

3. Absolute values of the housing vibration 
predicted by the global dynamic model did not 
always agree with measured values. 

4. The characteristics and trends of the housing 
vibration spectra predicted by the global dynamic 
model are the same or very similar to those found 
in the experimental data. 
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TABLE I. - TEST GEAR PARAMETERS 


Gear type Standard involute, full-depth tooth 


Number of teeth 28 

Module, mm (diametrial pitch in.* 1 ) 3.174(8) 

Face width, mm (in.) 6.35(0.25) 

Pressure angle, deg 20 

Theoretical contact ratio 1-64 

Driver modification amount, mm (in.) 0.023 (0.0009) 

Driven modification amount, mm (in.) 0.025(0.0010) 

Driver modification start, deg 24 

Driven modification start, deg 24 

Tooth-root radius, mm (in.) 1.35(0.053) 

Gear quality AGMA class 13 

Nominal (100 percent) torque, N-m(in.-lb) 71.77(635.25) 


TABLE 2.— COMPARISON OF 
EXPERIMENTAL MEASURED 
AND ANALYTICAL MODELED 
NATURAL FREQUENICES 


Experimental, 

Hz 

Analytical, 

Hz 

Difference, 

percent 

658 

658 

0 

1049 

1006 

-4.1 

1710 

1762 

3.0 

2000 

2051 

2.6 

2276 

2336 

2.6 

2536 

1536 

0 

2722 

2752 

1.1 

2962 

3012 

1.7 
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Figure 4.— Gearbox experimental mode shapes. 
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Major modes of housing 




Figure 5.— X-direction experimental and analytical vibration frequency spectra of 
the gearbox 
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Major modes of housing 




(b) Analytical. 


Figure 6. — Y-direction experimental and analytical vibration frequency spectra of 
the gearbox. 
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Major modes of housing 




7 00- Hz gear- \ Frequency, Hz 

mesh frequency - 1 


(b) Analytical. 

Figure 7.— Z-direction experimental and analytical vibration frequency spectra of 
the gearbox. 
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Mode: 1 

Frequency * 658 


Mode: 2 

Frequency * 1 006 



Mode: 3 Mode: 4 

Frequency « 1 762 Frequency « 2051 



Mode: 7 


Mode: 8 


Frequency * 2752 


Frequency » 301 2 


Figure 8,— Gearbox analytical mode shapes. 
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Displacement, X-direction, In. 

(a) Driving rotor orbit, (b) Driven rotor orbit. 

Figure 9. — Orbital motion of rotor during slow roll. 
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